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Abstract 

The energy of a graph is the sum of the moduli of the eigenvalues of its adjacency 
matrix. We study the energy of integral circulant graphs, also called gcd graphs, which 
can be characterized by their vertex count n and a set T> of divisors of n in such a 
way that they have vertex set Z n and edge set {{a, b} : a,b £ Z n , gcd(a — b,n) S V}. 
Using tools from convex optimization, we study the maximal energy among all integral 
circulant graphs of prime power order p s and varying divisor sets T>. Our main result 
states that this maximal energy approximately lies between s(p — l)^" 1 and twice 
this value. We construct suitable divisor sets for which the energy lies in this interval. 
We also characterize hyperenergetic integral circulant graphs of prime power order and 
exhibit an interesting topological property of their divisor sets. 
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1 Introduction 

Concerning the energies of integral circulant graphs, an interesting open problem is the char- 
acterization of those graphs having maximal energy among all integral circulant graphs with 
the same given number of vertices. The goal of this paper is to establish clarity concerning 
this question, for integral circulant graphs of prime power order, by showing how to construct 
such graphs with a prescribed number of vertices whose energy comes close to the desired 
maximum. In the course of this, we approximately determine the maximal energy itself. We 
rely on a closed formula for the energy of an integral circulant graph with prime power order 
that was established in |18|. 



1 



A circulant graph is a graph whose adjacency matrix (with respect to a suitable vertex 
indexing) can be constructed from its first row by a process of continued rotation of entries. 
An integral circulant graph is a circulant graph whose adjacency matrix has only integer 
eigenvalues. The integral circulant graphs belong to the class of Cayley graphs. By a result 
of So [20], they are in fact exactly the class of the so-called gcd graphs [20], a class that 
originally arose as a generalization of unitary Cayley graphs. The gcd graphs have first been 
described by Klotz and Sander in [TD] and further studied e.g. by Basic and Ilic [3J, 
[7]. The way the gcd graphs are defined serves us well, so throughout this paper we shall 
make use of this particular perspective of perceiving integral circulant graphs. Given an 
integer n and a set T> of positive divisors of n, the integral circulant graph ICG(n, T>) is 
defined as the corresponding gcd graph having vertex set Z n = {0, 1, . . . , n — 1} and edge set 
{{a, b} : a, b G Z n , gcd(a — b,n) G V}. We consider only loopfree gcd graphs, i.e. n V. 

The energy E(G) of a graph G on n vertices is defined as 

n 

E(G) = J2\M, 
t=i 

where Ai, . . . , A n are the eigenvalues of the adjacency matrix of G. This concept has been 
introduced several decades ago by Gutman [5], and with slight modification it can even be 
extended to arbitrary real rectangular matrices, cf . [13] and [9] . 

There exist many bounds for the graph energy, see Brualdi [1] for a short survey. One 
example is the bound 

71 

E(G)<-{V^+l) 

due to Koolen and Moulton [11] for any graph with n vertices. There exist infinitely 
many graphs that achieve this bound. If we consider only the class of circulant graphs, then 
the question arises how close one can get to this bound. Shparlinksi [19] has given an 
explicit construction of an infinite family of graphs that asymptotically achieves the bound. 

Another well-known result is due to Balakrishnan, who gives an upper bound B = 
k + \/k(n — l)(n — k) for the energy of a /c-regular graph on n vertices (see [2]). Li ET AL. 
[12] have shown that for every e > one can actually find infinitely many /c-regular graphs 
G such that > 1 — e. 

There has been some recent work on the energy of unitary Cayley graphs, which are 
exactly the gcd graphs with T> = {1}. Let us abbreviate £(ri,T>) = E(lCG(n,T>)) and let 
n = pi 1 ■ ■ -p s k k . Then, in the context of gcd graphs, the following result has been obtained 
by Ramaswamy and Veena [TJ] and, independently, by Ilic [7]: 

£(n,{l}) = 2V(n), 

where <p denotes Euler's totient function. 

Ilic [7J has slightly generalized this to some gcd graphs that are not unitary Cayley 
graphs: 

E(n, {l,pi}) = 2 k ~ 1 p i (p(n/p i ), provided that Sj = 1, 

E(n, {Pi,Pj}) = 2 fc y?(n), provided that si = . . . = Sk = 1. 
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In [18] the authors proved an explicit formula for the energy of ICG(p s , T>) for any prime 
power p s and any divisor set V = {p ai ,p a2 , . . . ,p ar } with < cii < a 2 < ■ ■ ■ < a r < s — 1, 
namely 

tr—l r \ 
p.-i r _(p_i)j^ j2 p s - a ^- i \. (i) 
^=1 i=k+i j 

The study of energies is usually linked to the search for hyperenergetic graphs. A graph G 
on n vertices is called hyperenergetic if its energy is greater than the energy of the complete 
graph on the same number of vertices, i.e. if E(G) > E(K n ) = 2(n — 1). Initially, the 
existence of hyperenergetic graphs had been doubted, but then more and more classes of 
hyperenergetic graphs were discovered. For example, Hou and Gutman show in [6] that if 
a graph G has more than 2n — 1 edges, then its line graph L(G) is necessarily hyperenergetic. 
Consequently, L(K n ) is hyperenergetic for all n > 5, a fact that seems to have been known 
before. 

Work by Stevanovic and Stankovic [2T] indicates that the class of circulant graphs 
contains a wealth of hyperenergetic graphs. Although integral circulant graphs are quite rare 
among circulant graphs (cf. [1]), the subclass of integral circulant graphs still exhibits many 
hyperenergetic members. For example, it has been shown by Ramaswamy and Veena 
[13] that almost all unitary Cayley graphs on n vertices are hyperenergetic. The necessary 
and sufficient condition is that n has at least 3 distinct prime divisors or that n is odd in 
case of only two prime divisors. Consequently, there exist no gcd graphs ICG(p s ,V) with 
T> = {1} that are hyperenergetic. However, for less trivial divisor sets it is also possible to 
find hyperenergetic gcd graphs on p s vertices. Some examples are given in [IB]. For instance, 
for p > 3 and s > 3, the choice T> = yields a hyperenergetic gcd graph. 

Not surprisingly, the class of graphs lCG(p s ,V) contains also non-hyperenergetic ele- 
ments, termed hypoenergetic. For the minimal energy £ m j n (n) of all integral circulant graphs 
with n vertices it has been shown in [18] that 

£min(p s ) = 2(p - l)/- 1 = E{K p s) - E{K p s-,). 

This follows directly from equation ([1]). The minimal energy is achieved exactly for the 
singleton divisor sets. 

The maximum energy of graphs ICG(p s ,P) is not as easily described. A classification of 
integral circulant graphs of prime power order p s with very small exponent having maximal 
energy has been provided in [18] , but it became clear that a general result as simple as in the 
case of minimal energy could not be expected. It will be the purpose of this article to clarify 
the structure of divisor sets imposing maximal energy on the corresponding gcd graph. Our 
main Theorem 14.21 states that the maximal energy among all integral circulant graphs of 
prime power order p s and varying divisor sets T> approximately lies between s(p — l)^ 8 ^ 1 
and twice this value. Tools from convex optimization will turn out to be the appropriate 
machinery to reach that goal. We shall compute bounds for the maximum energy and 
describe how to construct divisor sets for integral circulant graphs on p s vertices that have 
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near maximal energy. Along the way, we characterize hyperenergetic integral circulant graphs 
of prime power order and exhibit an interesting topological property of their divisor sets. 
Namely, the set containing all ordered exponent tuples corresponding to these divisor sets 
can be obtained by intersecting an integer grid with a suitable convex set. 

2 Preliminary definitions and results 

For any positive integer n, let 

E max (n) := max {£ (n, V) : V C {1 < d < n : d\ n}}. 

For given T> = {p ai ,p a2 , . . . , p flr } with < a\ < . . . < a r < s — 1, we have by ([1]), i.e. by 
Theorem 2.1 in [181, that 



£ (p s , V) = 2(p - ly- 1 (r - (p - 1)^(0!, . . . , Or)) , (2) 

where 

r— 1 r j 



for arbitrary real numbers Xi, . . . ,x r . In order to evaluate 'E max {p s ), our main task will be 
to determine 

Emaxty, r) := max {<£ (p s , V) : X> C {1 < d < n : d \ n}, |D| = r} 
as precisely as possible, given a fixed integer r. Therefore, we define for 1 < r < s + 1 

m p (s, r) := min {h p (a\, . . . , a r ) : < ai < a 2 < . . . < a r < s with a% G Z}. (3) 
It is then clear from ()2]) that 

£max(p s , r) = 2(p - l)p s_1 (r - (p - l)m p (s - 1, r)) . (4) 
Later on it remains to compute 

Emaxty ) = max {£ max (p s , r) : 1 < r < s}. (5) 
Proposition 2.1 Let p 6e a prime. Then 

(i) m p (s, 2) = ^ for all integers s > 1, and the minimum is attained only for a\ = and 
a 2 = s. 

(ii) m p (s, 3) = + ph + pS -\s/2] f or all integers s > 2. The minimum is only obtained for 
a\ = 0, a 2 = [s/2] (^or ; additionally, for a 2 = [s/2] + 1 if s is odd) and 03 = s. 

Proof. Proposition 3.1 in [18]. 



□ 

A set D C {1 < d < n : d | n} is called n-maximal if £(n, Z>) = £ m ax(^)- As a 
consequence of Proposition 12.11 and some other results in [18J, we obtained 



Theorem 2.1 Let p be a prime. Then 

(i) £max(p) = 2{p — 1) with the only p-maximal set T> = {1}. 

(ii) £ m ax(p 2 ) = 2{jp — l)(p+ 1) to^A £Ae only p 2 -maximal set T> = {l,p}. 

(Hi) £max(p 3 ) = 2(p — l)(2p 2 — p + 1) i/ie on/?/ p 3 -maximal set V = {l,p 2 }, except for 
the prime p = 2 for which T> = {1, 2, 4} z's a/so 2 3 -maximal. 

(iv) £ ma x(p 4 ) = 2(p — l)(2p 3 + 1) wi/i the only p 4 -maximal sets V = {l,p,p 3 } and V = 
{l,p 2 ,p 3 }. 



Proof. Theorem 3.2 in L , 

One can prove formulae for £ max (p s ) with arbitrary exponent s by using (j3J) and (J5J). As 
indicated in Proposition I2.1[ we need to choose integers < a± < a 2 < . . . < a r < s — 1 in 
such a way that 

i — 1 r 

/ lp (a 1) ..,o r ) = E E 

fe=i j=fe+i F 

becomes minimal. The choice of ai = and a r = s — 1 is clearly compulsory. 

The case r = 3 (cf. Prop. 12.1( h)) suggests to place ai, a 2 , . . . , a r _i, a r equidistant in the 
interval [0, s — 1]. A minor obstacle is the fact that the corresponding choice := i^ii— ij 
(1 < z < r) does not yield integral numbers as required. Taking nearest integers easily 
resolves this problem, but only at the cost of approximate instead of exact formulae. More 
seriously, it turns out that in general, even allowing real at, their equidistant positioning 
does not yield the desired minimum m p (s, r). The cases presented in Proposition 12. II do not 
yet exhibit this problem since it makes its debut for r = 4. An illuminating example can be 
found in the final section of [T8] . 

For the sake of being able to use analytic methods, we define for a prime p, a positive 
real number a and a positive integer r 

m p (a, r) := min {h p (a\, . . . , a r ) : < a.\ < a 2 < ■ ■ ■ < a r < a, a, € M}. (6) 

Observe that now the a« may be real numbers as opposed to integers in the definition of 
m p (s,r). It is obvious that for r > 2 

rh p (a, r) = min {h p (0, a 2 , ■ ■ ■ , a r _i, cr) : < a 2 < . . . < a r -i < cr, «i G R}- (7) 

Clearly, rh p (a, 2) = l/p CT , uniquely obtained for ot\ = 0, a 2 = a, and rh p (a, 3) = l/p°~+2/p a / 2 , 
uniquely obtained for a>i = 0, a 2 = a/2, = o (cf. Proposition 12 . ll fn) ) . 
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3 Tools from convex optimization 



In order to determine m p (a,r) in general it is crucial to observe that h p (xi, . . . , x r ) is a 
convex function. 

Proposition 3.1 Let r be a fixed positive integer, 6 ^ 1 a fixed positive real number and p 
a fixed prime. 

(i) The real function 

r r i ^ 

^'•••'^) : =EEn^ 

k=l i=k j=k 

is strictly convex on W . 

(ii) The function h p (xi, . . . , x r ) is convex on W . 
Proof. 

(i) Let (ui, . . . , u r ) 7^ (f i, . . . , v r ) be arbitrary real vectors. On setting 

i 1 * 1 

^M-n^ and ^ := n^-' 

j=k j=k 

we have Uk,k = b~ Uk and V^k = b~ Vk for 1 < k < r. Since < b ^ 1 and Uk ^ Vk for at least 
one k, we have Uk.k 7^ for that k. By the inequality between the weighted arithmetic 
and the weighted geometric mean, which is an immediate consequence of Jensen's inequality 
(cf. P], p. 1100, Thms. 17 and 18), we have U* ■ V 1 ' 1 < tU + (1 - t)V for all positive real 
numbers U and V and all < t < 1, with equality if and only if [/ = V. It follows that for 
all 1 < k < i < r 

uLv k y < tu k , + (i - 1^,, 

and strict inequality for at least one pair k, i. Hence 

r r 

g b {t(u 1 ,...,u r ) + (l-t)(v l ,...,v r )) = J]^^ 1 : 4 

k=l i=k 
r r r r 

k=l i=k k=l i=k 

= tg b (ui, . . . , u r ) + (1 - t)g b {vi, ...,v r ). 

This proves that g is strictly convex on W . 

(ii) Let (tii, . . . , w r ), (v i, . . . , f r ) G M r and < t < 1. Then by (i) 

h p (t(ui, ...,u r ) + (l- t)(vi, v r )) = 

= 9p(t{u 2 - Ui,« 3 - U 2 ,...,U r - Ur-l) + (1 - ~ Ul,«3 - «2, • • • , «r ~ «r-l)) 

< %(«2 - U h U 3 -U 2 ,...,U r - Ur-l) + (1 - *)fli,(u2 ~ V X ,V 3 - V 2 , ■ ■ ■ ,V r - V r -l) 

= th p (ui, ...,u r ) + (l- t)h p {v 1: ...,v r ), 
which shows the convexity of h p . 
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□ 



An easy corollary of ([I]) is a characterisation of the hyperenergetic gcd graphs of prime 
power order, namely that lCG(p s ,T>) with T> = {p ai ,p a2 , . . . ,p""~} and < a\ < a<i < . . . < 
a r < s — 1 is hyperenergetic if and only if 



(cf. Corollary 2.2 in [15]). As a first consequence of the convexity of h p we are able to 
refine this by showing that the set of hyperenergetic integral circulant graphs has a nice 
topological feature. Given a prime p and positive integers r < s, we define H(p s , r) as 
the set containing all (ai, . . . , a r ) G II with 0<ai<...<a r <s — 1 and the property 
that ICGfjr 5 , {p ai , . . . ,p ar }) is hyperenergetic. Then we can derive the following remarkable 
statement: 

Corollary 3.1 Let p be a prime and r < s positive integers. Then there is a convex set 
C C W such that H(p s ,r) = Cf]Z r . 

Proof. For fixed p, s and r, we define 



Since h p is convex on W by Proposition I3.1[ the so-called level set 

L := {(xi, . . . , x r ) E W : h p (xi, . . . , x r ) < c(p, s, r)} 

is also convex (cf. [17] . p. 8 and Prop. 2.7). Since 

K := {(xi, . . . , x r ) G W : < x x < x 2 < ■ ■ ■ < x r < s - 1} 

is obviously convex as well, the intersection C :— Lfl K has the same property. By (jSJ) we 
know that some (a 1; . . . , a r ) G Z' r lies in H(p s , r) if and only if < a\ < . . . < a r < s — 1 
and h p (ai, . . . , a r ) < c(p, s, r), hence "H(p s , r) = C fl Z r . ^ 

We shall use some further standard results from convex optimization. 
Proposition 3.2 Let f : U — >■ R be a strictly convex function defined on a convex set 



(i) If U is an open set then each extremal point of f is a minimum. 

(ii) If f has a minimal point on U then it is unique. 

Proof. The proofs of the assertions can be found in [16J, pp. 123-124, Theorems A 






U C W. 



and C, in [8], or in [T7], Thm. 2.6. 
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□ 



Our main tool for the computation of rh p (a,r) will be 



Proposition 3.3 Let r > 1 be a fixed integer. We define the real function 



r 



r 



f(x u ...,x r ) := x i 



k=l i=k j=k 



for (x u ...,x r )e [0, l] r . Let < p < 2~ r . Then 

min{/(xi, . . . ,x r ) : (xi,...,x r ) E [0, l] r , Xl ■ x 2 ■ ■ ■ x r = p} = (r + n(p,r)) • n(p,r), 

where p(p,r) :— v(p, r)/(l— v(p, r)) and x = v(p,r) is the unique real solution of the equation 
x r = p(l — x) 2 on the interval [0, 1]. The minimum obviously equals p for r = 1, and it is 
p + 2^pforr = 2. 

There is a unique minimizer for each r, namely 



PROOF. For r < 2, we have to deal with nothing more than quadratic equations, and 
in these cases all assertions follow easily from standard analysis. 

For r > 3, we use the method of Lagrange multipliers to obtain necessary conditions for 
local minima of f(x±, . . . , x r ) subject to the constraint x\ ■ x 2 ■ • • x r — p. Accordingly, let 



A necessary condition for a local minimum is that all partial derivatives F Xt := (1 < t < r) 
as well as F\ := §^ vanish at that point. We have for 1 < t < r 




In the special case r = 2, we have explicitly p(p, 2) = ^fp. 



F(x u ...,x r ,X) := f(x u ...,x r ) + X(p-x 1 -x 2 ---x r ). 






Hence 





k=l i=t j=k 



r 
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and we want to find all solutions (xi, . . . , x r ) of the following system of equations: 

( t r i 

JEEII^ = V (l<f<r), 

k=l i=t j=k V*) 



1 



X\ ■ • ■ x r — p . 

From now on we consider r and p to be fixed and abbreviate p := p(p, r) and v := v(p, r). 

Claim : One solution of fl9]) is given by X\ = x r = p and X2 = ■ ■ ■ = x r _i = u, where we have 
< v < p < 1. 

The real function h(x) := x r — p(l — x) 2 is strictly increasing on the interval [0, 1] with 
h(0) = —p and h(l) = 1. Hence x r = p(l — x) 2 has a unique solution on (0, 1), which is 
denoted by v. Since u r < p < 2~ r , we even know v < 1/2. This implies v < p, < 1. 

In order to show that (p, u, is, . . . , v, p) G (0, l) r satisfies (ED, we separate terms containing 
x\ or x r from the others in the double sum of and obtain for x\ — x r — p and X2 — ■ ■ ■ — 

X r -\ = v 



t r i r— 1 i t r t r—1 i 

k=l i=t j=k i=t j=l k=2 j=k k=2 i=t j=k 

r-1 t t r-1 



1 I 



i=t k=2 k=2 i=t 



v i-k+l 



Jjt-1 j/""" 1 t 1 * yt—k+1 yT— fc+1 



v — v v — V X — ^ 

p + p — + p — + > 



1 — V ' 1 — V ^— ' 1 — V 

k=2 

J ~ l 2u r - l + u r '- t 1 (v-v 1 v'- t+l 



V — LV -\- V 1 I V — V I/' lT - V 

= P + P ; + ; 1 ; 

1 — V 1 — V \ 1 — V \ — V 

= P + M 1 (M 1 - - 2 ^ _1 + O + - - ^ _t+1 + ^)) • 

(1 — vy 

Since p{l — v) = u, we conclude for . . . , x r ) = (p, v, . . . , 1/, p) that 

u-{l-u r - 1 ) p 



1 - u) 2 " = i>-i 



k=l i=t j=k 



Setting A := l/u r 1 , this reveals that (p, v, . . . , u,p) satisfies all the upper equations in 
The observation that 

2 r-2 V 

X\ ■ ■ • x r = p -v = — = p 

(1 - vf 

completes the proof of the claim. 

We now want to show that f(p, u, . . . ,u, p) is in fact a minimum subject to the constraint 
x\ ■ X2 ■ ■ • x r = p, and we shall see as well that (p, u, . . . , u, p) is the unique minimizer. By 
Proposition 13.1( 1) the function #2(2/1, • • • , Vr) = /(2~ yi , . . . , 2~ Vr ) is strictly convex for all 
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(yx, . . . , y r ) G W. Our claim has shown that (p, v, . . . , is, p) is an extremal point of / on the 
set {(xi, . . . ,x r ) G [0, l] r : X\- ■ -x r = p}. Therefore (//, is', ... , is',p') with 

p log 2 log 2 

is an extremal point of g% on the set 

U := . . . , y r ) G M r > : 2"* ■ ■ ■ 2"* = p} = . . . , y r ) G M r > : ^ + . . . + y r = a}, 

where 

a:= MlM> r . 

log 2 " 

The set U apparently is the simplex with vertices (a, 0, . . . , 0), (0, a, 0, . . . , 0), . . . , (0, . . . , 0, a), 
and therefore a convex subset of E> . It is immediately clear that (p', is', ... , is', p') does not 
lie on the boundary of the simplex U, in other words: the point belongs to the set U° of 
inner points of U. Altogether the function #2 is strictly convex on the open convex set 
U° and (p 1 , is', ... , is',p') is an extremal point of gi in U°. By Proposition I3.2H ) this point 
(p 1 , is', ... , is', p') is a minimal point of g 2 , and by Proposition 13.2( h) it is unique. Since log 2 
is strictly monotonic, the point (p,is, . . . , is, p) is the unique minimizer with respect to / on 
{{xx, . . . , x r ) G [0, l] r : x x ■ ■ ■ x r = p}. 

It remains to calculate the minimum. We obtain similarly as before 



r— 1 i r r r— 1 r— 1 i 

f(x 1 ,...,x r ) = Xl x r + j2H x+ J2Il x+ J2J2Il 

i=l j=l k=2 j=k k=2 i=k j=k 



Xj. 



By evaluating the geometric sums and using the identity p = is/(l — is), it follows that 

1 — z/' -1 If is 2 — is r * 

f(p, is,...,is,p) = p 2 is r ~ 2 + 2p ■ — + (r - 2)is - 



1 — V 1 — V \ 1 — V 

2u(l - is*- 1 ) is . v 2 - v T 

+ -4: ^ + z r - 2) 



(1 — is) 2 (I -is) 2 I -is " ' (l-isf 

2v-v 2 v 



+ Z {r - 2) 



{I -is) 2 I -is 

2v v 2 v 

+ TT-. ^ + ■ (r - 2 ) 



I -is ' {I -is) 2 ' I -is 
2p + p 2 + p(r — 2) = (r + p) ■ p. 



□ 



Corollary 3.2 Let p be a prime and r > 2 an integer. For a given real number 
a > (r — 1) log 2/ logp let x = v v {p, r) be the unique real solution of the equation p cr x r ~ 1 = 
(1 — x) 2 on the interval [0,1], and p p (a,r) := v v {a, r)/(l — v v {a,r)). Then 



(or, r ) = ( r - 1 + p p (cr, r)) ■ p p (cr, r). 
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and this value is exclusively attained by h p (at\, . . . , a r ) for ctj = <x,(cr, r) (1 < j < r) defined 
as ai(a,r) := 0, a r (a,r) := a and 

log (u p (a, r)' 1 ) log (z>„(cr, r) _1 ) 

aj ((T,r):= &V 7 V J +(3-2) & K M ; r > 3; 2 < j < r - 1 . 10 
logp logp 

Proof. Let < ai < a 2 < . . . < a r < ff be arbitrary, and set yj := atj+i — atj for 
1 < j < r — 1. Hence a-j — = yt + yu+i + • • • + y%-\ for 1 < k < i < r. This implies 



r— 1 r | 

fc=l i=fc+l r 

r— 1 r «— 1 . r— 1 r— 1 j ^ 

fc=l i=k+l j=k k=l i=k j=k 

On setting := for 1 < j < r - 1, we have h p (ai, . . . , a r ) = /(^i, . . . , x r _i) for the 
function / as defined in Proposition 13.31 By hypothesis, r — 1 > 1 and (x\, . . . , x r _i) G 
[0, l] r_1 . Now we search for conditions to be imposed on the ay in order to hit the minimum 
rh p (a,r). First of all, we necessarily have a± = and a r = a according to (171). Hence 

11 1 

X\ ' X2 ' ■ ■ ■ ' X, — 1 — ; ; — — . 

pVl+.-.+Vr-l pOtr—O-l 

Again by hypothesis 

1 1 
< p := — < -. 

Applying Proposition I3.3j we conclude that 

m p (a, r) = (r - 1 + fi(p, r - 1)) • fi(p, r - 1) = (r - 1 + Jx p {a, r)) ■ p, p (a, r), 

where this minimum is exclusively obtained for = fi p (a, r) and x 2 = x 3 = . . . = 

x r -2 = P p (a,r). This yields for the given values of the Xj 

1 1 1 ~ / \ 

— = = — = x 1 = u p {a, r), 

pa 2 pct 2 -a 1 pyx ± r~p\ i n 

hence a 2 = log(l//2 p (cr, r))/logp, and for 2 < j < r — 2 

1 1 . , , 



pay+i-ay pVj 

which implies (fTUj) . 



xj = v p {a,r), 



□ 
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4 Integral circulant graphs with maximal energy 

Up to this point, all we have done with respect to general integral circulant graphs with 
maximal energy refers to real parameters atj in h p (ai, . . . , a r ). As a consequence, we have 
the following upper bound for £ max (p s , r), but we are left with the task to find out how close 
we can get to the "real maximum" if we restrict ourselves to integral parameters 
as required by our problem. 

Theorem 4.1 For a prime p and integers 2 < r < s, we have 

£max(p s , r) < 2(p - ly- 1 (r - [p - 1) (r - 1 + jx p (s - 1, r))fi p {s - 1, r)) , 

where jl p is defined in Corollary VJ. 6 ^ 

Proof. By Corollary 13. 2l and the definitions of m p (s— 1, r) and m p (<7, r), we immediately 
have for any integer s > r — 1 

m p (s - 1, r) > (r - 1 + fi p (s - 1, r)) ■ fi p (s - 1, r). 

Now our theorem follows at once from this and (jl]). 

□ 

The first step we take towards integrality of the parameters is to approximate the numbers 
p, p (s — l,r), u p (s — l,r) and the corresponding ctj(s — l,r), all defined in Corollary 13.21 by 
simpler terms. 

3-1 

Proposition 4.1 For a prime p and integers 3 < r < s, let 5 := p r - 1 . Then we have 

(i) 6<fi p {s-l,r)<6 + £- s ; 

(ii) 8-&<v p (8-l,r)<5<±; 

(m) < log (z> p (s - 1, r)- 1 ) - g \ogp < 

(iv) < log (p, p {s - 1, r)" 1 ) - logp < 0; 

(v) \aj(s - l,r) - (j - 1)^| < ^ /or 1 < j < r. 

PROOF. (i) It follows from the definition of jl := fl p (s — l,r) in Corollary 13.21 that it 
satisfies the identity p s ~ 1 /i' r ~ 1 = (1 + /l) r ~ 3 , clearly implying < jl < 1. For r = 3, this 
means that jl = 5. For r > 4, we obtain by virtue of binomial expansion 

ji = 5{l + ji) r — =S + 5j2[ r k 1 

k=i ^ 
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where the infinite series has alternating decreasing terms. Hence 

y ^ 

< ft — 5 < 5 ■ ■ ft < 5 ft, 

r — 1 

and consequently ft < 5/(1 — 5), which implies (i). 

(ii) Since the real function x i->- x/(l +x) is strictly increasing for x > 0, we obtain by (i) 

A ^ 6 < » , A 

1 + 5 1 + 5 - l + 1 + 5+^ 

The definition in Corollary 13.21 yields that v := P p (s — l,r) = ft /{I + ft), which proves our 
claim. 

(iii) By (ii), we have < v < -. Taking logarithms in the identity p s ~ x v r ~ x = (1 — ii) 2 
(cf. Cor. 13. 2p . we obtain 

1 s — 1 2 

log- = -logp -log(l-z/). (11) 

v r — 1 r — 1 

Since v < 1/p, the Taylor expansion of log(l — v) yields 

k=l k=0 

Inserting this into (TTTj) . we get 



; — i z — n \ / ^ 



i 1 s-1 3 
< log — log p < 



v r — 1 (r — l)p 

(iv) For the numbers ctj := ay(s — 1, r), as defined for 1 < j < r in Corollary 13. 2[ we have 

{ log(V^) f or „• _ ^ anc [ j = r _ i 
« for2< j<r -2. <12) 

This is trivial except for j = r — 1, where it follows from the identities ft = vj{\ — v) and 
ps-i~r-i _ ^ _ therefore, 

^ 21og(l//2) (r-3)log(l/z>) 

s — 1 = a r — ai = > (a 7+ i - a,) = — 1 : , 

^— < logp logo 

i=i 

hence 

log logp = — — logp - log - . 

ft r — 1 2 \r — 1 z// 

Combining this with the bounds found in (iii) completes the argument. 



v 



) By the definition of the ctj, we obtain for 2 < j < r — 1 



.s-1 /log \ s-l\ , , /logi s-l\ 

^ s — 1, r) — (j — 1 = -j — + j - 2 

r — 1 \ logp r — 1 I \iogp r — 1 J 
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From (iii) and (iv) it follows that 



3 , .si, 3 3 

< aj {s - i, r ) - - 1)— r < - 2 ) ?I ^„i„,„ < 



2p log p r — 1 (r — 1 )p log p p log p 

which implies (v) in these cases. Since a% = and a r = s — 1, the inequality is valid for all 
3- 

□ 

Proposition 14. lf v) reveals that picking the atj for j = 1, . . . , r well-spaced in the interval 
[0, s — 1], i.e. ctj := (j — 1)^Et (see concluding remarks of section 2), is close to best possible. 
Since it is our task to find integral aj in optimal position, it suggests itself to choose the 
cij as nearest integers to the a>j(s — l,r) (as defined in Corollary 13. 2p or to the numbers 
(j — l) 2 ^) which does not make much of a difference by Proposition 14. lf v). Anyway, we 
shall take aj = \\atj\\ (1 < j < r) with the nearest integer function || ■ || and have to accept 
variations between aj and aj in the range from — | to \. Finally, we shall try to maximize 
the energy with respect to r. 

We now show that our integral minimum m p (s — 1, r), as defined in ([3]), can be bounded 
by the real minimum m p (s — 1, r), introduced in (Q. In general, that is to say in worst cases, 
we cannot expect to lose less than a factor p between the two minima, taking into account 
that the shifts from real numbers aj to integral parameters aj, varying over an interval of 
length up to 1, have to be executed in h p , i.e. in the exponent of p. 

Proposition 4.2 Let 3 < r < s be given integers, and let p be a prime. 

(i) Let («!,..., a r ) G W be the unique minimizer of h p determined in Cor. VS. ''A Then 
(ai, . . . , a r ) G Z r with the nearest integers aj := \\aj\\ (1 < j <r) has the property 



hJai, ...,a r ) < 



4 ■ m2(s — 1, r) for p = 2, 
p-m p (s — l,r) forp>3. 



(ii) We have 7712(5 — 1, r) < iri2(s — 1, r) < 4 ■ m2(s — 1, r). 

(iii) For any prime p > 3, we have 

m p (s — 1, r) < m p (s — 1, r) < p ■ rh p (s — 1, r). 

Proof. The lower bounds in (ii) and (iii) are trivial, and the upper bounds follow 
immediately from the definition of m p (s — l,r) in (j3]). Hence it suffices to prove (i). By 
Corollary 13.21 we have aj = <x,(s — 1, r) for 1 < j < r and 

h p (ai, ...,a r )= rh p (s - 1, r). 
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Since v := v p (s — l,r) < 1/p by Proposition I4.1f ii). we have \og(l/v) > logp. From 
Proposition l4.1( i) it follows for p > 3 that jl := jX p (s — l,r) < 1/a/p, hence log(l//j) > |logp. 
By use of (fT21) . these inequalities imply that a 2 > 5, ol t -\ < s — 1 — \ an d > ctj + 1 for 
2 < j < r — 2. Moreover, a\ = and a r = s — 1. Therefore, the nearest integers aj := ||<x,||, 
1 < j < r, are pairwise distinct, forming a strictly increasing sequence. We have aj = aj + 5j 
for suitable real numbers 5j satisfying \Sj\ < 1/2 (1 < j < r) and obtain 



1 

/i p ( ai ,...,a r ) = ^ ~Z 

k=l i=k+l y 

r—lr r—1 r 

= V V - = V V — — 

Z_-^ Z-^l r)(oii—Si) — (a k —5 k ) Z-^/ Z—J ff&k—Si rPi~ a k 

k=l i=k+l y k=l i=k+l J 

r—1 r ^ 



k=l i=k+l y 

For the prime p = 2 the above proof has to be modified, since possibly a 2 < |. In this case 
we choose a 2 = 1 and a r _i = s — 2. As before, ay+i > ctj + 1 for 2 < j < r — 2. Hence we 
can select each aj, 3 < j < r — 2, as one of the neighbouring integers of <x,- in such a way 
that ax < &1 < . . . < a r . It follows in this case that aj = aj + 8j for suitable real numbers 
5j satisfying \Sj\ < 1 (1 < j < r). Consequently 

r—1 r 1 1 r ~ * r X 

h 2 (a u ...,a r ) = J2 E o^o^^E E ^p^ = 4 -^(3-l,r). 



2 <5 fc- (5 i 2 ai ~ ak ~ ^—^ 2 C 
fc=i t=jfe+i fc=i t=jfc+i 



□ 



Remark. The factor p (or 4 in case p = 2, respectively) we lose between the real minimum 
777.2(5 — l,r) and the integral minimum m p (s — l,r) according to (ii) and (iii) reflects the 
hypothetical worst case scenario where each aj differs from aj by |. In practice, the factor 
between the two minima will be substantially smaller in almost all cases. 



In Theorem 12.11 the maximal energy £ max (p s ) as well as the corresponding p s -maximal 
sets are given for all primes p and each s < 4, and could be determined quite easily for other 
small values of s by (J2J), i.e. Theorem 2.1 in [TB]. The inequality 

is an immediate consequence of Theorem 14.11 The following result shows that this trivial 
upper bound lies close to the true value of E*^^ 8 ). 

More precisely, part (ii) of the following Theorem 14.21 provides the explicit construction 
of a divisor set T>o such that the energy of the graph lCG(p s , T>q) falls short of the maximal 
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energy T, max (p s ) among all integral circulant graphs of order p s essentially by a factor less 
than 2. The remark preceding Proposition 14.21 explains why we cannot expect to find a more 
precise lower bound in general. However, the reader should be aware of the fact that we lose 
a much smaller factor than 2 between upper and lower bound for £max(p s ) m m ost cases (cf. 
the remark following Prop. 14. 2p . We shall comment on this at the end of the section. 
Bound by the tradition of number theory, log will denote the natural logarithm. 

Theorem 4.2 Let p be a prime and let s be a positive integer, 
(i) We have 

log logp 



Q.(s-l) 1 



logp 



<zL x (p s )<c-(s-i) l- 



log logp 
logp 



+ 1, 



(13) 



where C — 1 for all p > 3 and C_ = \ for all p > 17 as well as for 3 < p < 13 in 
case s < 6. Only for small values of p, we have exceptional constants C = C(p) and 
C = C(p), namely C{2) = 0.328, C(2) = 0.118, and in case s>7 



C{p) 



0.030 ifp = 3, 

0.233 ifp = 5, 

0.337 ifp = 7, 

0.442 if p= 11, 

0.473 ifp=U. 



(ii) Let r be the integer uniquely determined by 

s 



1 s 

<r < 



D(p) 



where 



Dip) :-- 



4.09184 

2(1 + l2£]2S£^ 



+ 1, 

forp = 2, 
for p > 3, 



and define T>o = {p» a ^ s 1 ' r °)W : j = 1 ; . . . , ro}. For p — 2, s > 11 and for p > 3, s > 7, 
t/ie energy of the graph lCG(p s ,T>q) lies in the same interval as the one established for 
£ m axG» s ) m ([0 



Proof. By Theorem 12.11 we have for all primes p 

1 



2 + 4, 



for s 
for s 
for s 
for s 



(14) 



4. 



We leave it to the reader to check that each of these values lies within the respective bounds 
stated in (TLB"]) . We may therefore assume s > 5 in the sequel. 
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We shall first prove the upper bound in (|T3i) . By virtue of ([5]), it suffices to show that 
£max( ^ r) <C.(s-l) (l- l °f l0gP )+l. (15) 



(16) 



2(p — l)p s 1 \ logp 

is satisfied for all 1 < r < s. We distinguish three cases. 
Case Ul : r — 1. 

By Corollary 2.1 (i) in [18], we have £ m£K (p s , 1) = 2(p — l)p s_1 for all p, which implies f[To"j) 

immediately. 

Case U2 : r = 2. 

It follows from (J4]) and Proposition I2.1( i) that 

£max(p S , 2) = 2(p - I)/" 1 U-ip- 

for all p. Since s > 5, our upper bound in (1X5]) is valid in this case. 
Case U3 : 3 < r < s. 

By Theorem 14.11 and Proposition I4.1f i). we have 

£ max (pV) < 2(p- ljp- 1 (r- (P- 1 ^- 1 ) ) (i 7 ) 

V p^ 1 / 

for all p. Therefore, we study for fixed p and s the real function 

(p-lXs-l) 
#(x) := x ^ 

on the interval 3 < x < s with boundary values 

g(3) = 3-^=^ and «j( a ) = — + 1 . (18) 
p— p 

For a maximum of # at Xq, say, with 3 < xq < s the derivative 
vanishes necessarily. Substituting y := hence y > 1, we obtain the condition 

p2/0 

l + y logp= r (19) 

p — 1 

for y := -^r. Since x = 1 + (s — l)/yo, we conclude for 3 < x < s 

9(*) < 9(x ) = — ( 1 - ^) + 1. (20) 

Case U3.1 : p = 2. 
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For p — 2, equation (TT9]) has no solution > 1> i-e. in that case the maximum of (7 is 
attained at yo = 1, that is for xo = s > 5, which follows by comparison of the boundary 
values in (lisp . Therefore we have in case p = 2 

s — 1 s — 1 
#(z) < S-(s) = s ^— = -y- + 1. 

By (lTTj) . this immediately implies ( TT5"j) . 
Case U3.2 : p > 3. 

Now (Tl9l) has a unique solution ?/o in the interval 1 < yo < 2, corresponding to the unique 
maximum of g. By a few steps of Newton interpolation we obtain for instance that yo ~ 1.527 
for p = 3 and ?/o ~ 1-673 for p = 5. We shall verify that 

9W <( S -l)(l-te)+l (21) 

on the interval 3 < x < s. This follows easily for p = 3 and p = 5 by inserting the respective 
values of ?/o given above into (|20|) . For each other fixed prime p > 7, we define the real 
function 

w{y) = w p {y) := - y \ogp - 1 

p — 1 

for all y > 1. By (TIP]) we know that y > 1 satisfies w(?/o) = 0. Since the derivative 

w'(y) = (^yZi ~ X ) logp 

is positive for y > 1, the function is strictly increasing. For 

logp 1 
^ p log p — log log p log p ' 

which is greater than 1 for p > 7, we have 

/ \ 1 , — Jsp, — -1 (logp) 2 
w{y p ) = p log p-iog logp . e 



p — 1 log p — log log p 



1 1 I log logp (lOgp) 

-p log p-log logp 



2 



e(p — 1) logp — log logp 

P log logp 

< — log p-log logp _ logp < , 

e(p-l) 

where the final inequality is shown to hold for all primes p > 7 by simply taking logarithms 
in 

P log logp 

-p log p-log logp < lOgp. 



e(p- 1) 

Since w(y) is strictly increasing on y > 1 and w(y p ) < 0, but w(yo) = 0, it follows that 
y p < yo. By definition of y p , this inequality implies 

logp < ( 1 - loglogp J (1 + y logp). 
V. logp J 
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Multiplying with yo and dividing by (1 + yologp), we obtain by (1T91 

p-1 1 / loglogp 

1 ZT = 1 - i < 1/0 1 ~ 



p yo l + ?/ologp V logp 

Inserting this into (|20|) . we have verified (l2"Tj) . By ( 1T7|) . the proof of ffTBT) is complete. Hence 
the upper bound in (fTBT) holds in all cases. 

Now we turn our attention to the lower bound for £m ax G° s ) an distinguish several cases 
and subcases. 
Case LI : p > 3. 
Case Ll.l : s < 4. 

The lower bound in (TTBj) has already been verified for s < 4 in (|14p . 
Case LI. 2: s = 5. 



Picking r = 2, we use (1X6]) once more and obtain 

^<P 5 )>?=^ = 2-( P -l)A>2fl '° gl0gP 



2(p — l)p 4 p 4 \ \ogp 

for all p > 3, which proves the lower bound of ffTBT) in this case. 
Case LI. 3 : s = 6. 

It follows from (J3J) and Proposition I2.1( ii) that 

* / fis £max(p 6 ,3) . . / 1 1 1 \ 5 / log log V 



2(p — l)p 5 ' \p 2 p 5 p 3 J 2 \ logp 

for all p > 3, and again the lower bound in ({TBI is confirmed. 
Case L1.4 : s > 7. 

We choose the integer r according to the inequality 

s — 1 s — 1 , „ 

^- < r < — - + 1, (22) 

ZiJ.Jp £±jp 

where 

log log p 

Li v :— l H . 

logp 

A simple calculation reveals that for s > 7 and all primes p > 3 the expression on the 



left-hand side of (122]) is always greater than 2. Consequently, 3 < ro < s — 1. By Corollary 
13. 2[ we have 

m P - 1, r ) = (r - 1 + /2 p (s - 1, r )) ■ ft p (s - 1, r ), (23) 

and from Proposition I4.1( i) and the definition of r , we get p, p (s — 1, r ) < 5 + <5 2 / (1 — <5) for 
5 = p-(«-i)/(»-o-i) < p -2L P , Hence 

/ip(s - l,r ) < 
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By (123]) and (122]) we now have 

s-1 1 



m p (s-l,r ) < 



+ 



1 + 



1 



< 



2L p p 2L p \ p 2L p — 1 

s — 1 1 \ » 
+ 



p2Lp 



l + 



,2 



because for p > 3 



1 + 



p2L p — Y 



2p 2L pL p p iL p J p — 1 



< 1 + 



p2L p — ]^ 



< 1 



< 



p 2 (\ogp) 2 — 1 

p 



(p — l)(p + 1) J p-1 
Since p > 3, we obtain by Proposition I4.2( iii) that 



,(s-l,r ) < p-m p (s-l,r ) 



< 



s - 1 



2p 2L pL p p 4L p J p — 1 

s-1 1 
+ 



2(logp) 2 L p p 2 (logp) 4 L p y p-1' 
According to (gD and ([22]), it follows that 

£max(p S ,r ) 



2(p - l)^- 1 



> 



> 



r -(p- l)m p (s - l,r ) 

s - I / s - 1 1 

+ 



2L p \2(\ogp) 2 L p p 2 (\og P Y 
s - 1 

2Lr, 



1 



It is easy to check that for all primes p > 17 

(log log p) 2 > 1 

and that for 3 < p < 13 

(log log p) 2 > 1 - 



(logp) 2 / p 2 (logp) 4 



3(p logp) 2 



1 



3(p\ogp) 2 



c p {log p) 2 L 



with constants c„ defined by the following table: 



p 




3 


0.859 


5 


0.375 


7 


0.214 


11 


0.073 


13 


0.033 
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Setting c p := for all primes p > 17, the fact that s > 7 implies for all p > 3 

L 2L 
(loglogp) 2 + Cp\ogpL p > 1 + p > 1 + 



3(plogp) 2 (s — l)(p\ogp) 2 

Dividing by (logp) 2 , adding 1 on both sides and rearranging terms yields 



1 1 2L P ^ log log 2 CpLp 



(logp) 2 (s — l)p 2 (logp) 4 \ logp / logp 

log log p\ T c p L p 



logp / logp 
Dividing by L p and multiplying with (s — l)/2 implies 

s - 1 f 1 \ 2 s - 1 / log log p + c p 



2L p \ (logp) 2 / (s — l)p 2 (logp) 4 2 \ logp 
Inserting this inequality into (I24j) yields 



(p s , r ) s - 1 / log log p + c p 



2(p - l)^" 1 ' 2 V. logp 



which completes the proof of the lower bound in ( Tl3[) for p > 3. At the same time, our 
construction combined with Proposition 14. 2 ( i) reveals the truth of statement (ii) for primes 
p > 3 and s > 7. 

Case L2 : p = 2. 
Case L2.1 : s < 4. 

The lower bound in ffTB"]) follows from (jl4j) . 
Case L2.2 : 5 < s < 10. 

Picking r = 3, it follows from (T4]) and Proposition I2.1f ii) that 

rr* (QS\ > £ max(2 S , 3) / 1 1 1 



S-ll 



It is easy to check that the last term becomes minimal for s = 5. Hence 

£; ax (2 s )>3-^>^( S -i) 

for all s in the given range. This confirms the lower bound in (fT3"l) for these values of s. 
Case L2.3 : s > 11. 

Let c\ > 1 be the unique real number satisfying c 2 — 6ci + 5 = 4cilogci, i.e. 17.0517 < 
ci < 17.0518, and let c 2 := ^gf , thus 4.09184 < c 2 < 4.09186. We choose the integer r 2 
according to the inequality 

s — 1 s — 1 , . 
< r 2 < + 1. 25 

c 2 c 2 
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Apparently, the expression on the left-hand side of (1251) is always greater than 2 for s > 11. 
Consequently, 3 < r 2 < s — 1 . By Corollary I3.2[ we have 

m 2 (s - 1, r 2 ) = (r 2 - 1 + fi p (s - 1, r 2 )) ■ /2 p (s - 1, r 2 ), (26) 

and from Proposition I4.1f i) and the definition of r 2 , we get 

1 1 If 1 \ 

By (T26D and ([25]) we obtain 

m 2 (s-l,r 2 )< (— + ^(l + J-r)V^[l 



c 2 2 C 2 V 2C2 - 1 / / 2 C 2 V 2C2 - 1 / ' 

Proposition I4.2( ii) implies that m 2 (s — l,r 2 ) < 4 • m 2 (s — l,r 2 ), and with (T4]) and (|2"5]) we 

get 

£ max(2 ) > ^ = r 2 -m 2 (s - l,r 2 ) 

s- 1 /s - 1 1 / 1 \\ 1 / 1 
> 4 + — - 1 + ^- ■— 1 + 



c 2 V c 2 2 C 2 V 2 C 2 - 1 / J 2« V 2 C 2 - 1 



c 2 V 2 c 2 V 2 c 2 - iy y 2 2c 2 V 2 c 2 - 1 

> IZi- 0.01553 >2(.-l) 

for s > 11. This proves the lower bound of (Tl3i) and completes part (i) of Theorem 14.21 This 
time our construction combined with Proposition 14. 2f i) shows (ii) for p = 2 and s > 11. 

□ 



As an example, the following table illustrates the previous theorem for s = 17 and several 
small values of p. Note that r will eventually become 8, roughly for p > 10 10 ). 



n 


^0 


V 


lower 


^max( n ) 


upper 


3 17 


8 


(0,2,5,7,9,11,14,16) 


0.439 


6.652 


15.630 


5 17 


7 


(0,3,5,8,11,13,16) 


2.626 


6.547 


12.269 


? 17 


6 


(0,3,6,10,13,16) 


3.547 


5.927 


11.526 


ll 17 


6 


(0,3,6,10,13,16) 


4.493 


5.969 


11.164 


13 17 


6 


(0,3,6,10,13,16) 


4.789 


5.978 


11.124 


17 17 


6 


(0,3,6,10,13,16) 


5.059 


5.987 


11.119 


23 17 


6 


(0,3,6,10,13,16) 


5.084 


5.993 


11.169 



Remarks. 
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(i) The proof of Theorem 14. 21 shows that we lose a factor 2 between lower and upper bound 
for Emax(P s ) m case P — 17 (and similarly for smaller p) mainly due to the fact that we 
lose a factor p between fh p (s — 1, r) and m p (s — 1, r), which however happens only as an 
extremely rare worst case event (cf. Proposition 14.21 and the preceding and subsequent 
remarks). A staightforward adaptation of the method introduced in the proof for the 
lower bound of (TIB"]) implies the following: 

Assume that for some fixed sufficiently large p and s we have m p (s — 1, r) < p 1 ■ m p (s — 
l,r) with some positive 7 < 1. Taking 

8 + 1 8 + 1 , 1 
7 — < r < : h 1 

instead of ( 125|) . we obtain 

1 + 7 V. logp / V. logp J 

This reveals that, the smaller the difference between integral and real maximum is, 
the better our bounds are. In the extreme case where m p (s — l,r) = m p (s — l,r), 
i.e. 7 = 0, lower and upper bound differ only by 1, and the gcd graph with the 
corresponding divisor set T>q has maximal energy, because the energy is an integral 
number, but the two bounds are not. 

(ii) It should be noted that even the lower bound of ( TTBl already implies hyperenergeticity 
in most cases. A straightforward calculation shows this for e.g. p > 3 and s > 5. 



5 Conclusion and open problems 

Given a fixed prime power p s , we have provided a method to construct a divisor set T>q 
with the property that t E(p s ,V ) > ^ < E raax (p s ). In most cases we expect < E(p s ,V ) to lie 
much closer to £ ma x(p s ) than our worst case inequality guarantees. But since we have used 
the "real" maximum for reference it may not be expected to get hold of the "integral" 
maximum in general, using an analytic approach. The convexity properties of the function 
h p also suggest that a divisor set T> max with £(p s ,P max ) = £ max (p s ) can be found in the 
"neighborhood" of our T) . Given an explicit integer p s it should not be too difficult to 
determine £ ma x(p s ) precisely by comparison of a few candidates for a P max "near" Do- 
Let us conclude this section by posing the challenge of finding similarly accessible bounds 
on £ max (n) for integers n which have different prime factors. Even for n = pj 1 ^ 2 with primes 
Pi 7^ P2 and arbitrary divisor sets a closed formula for the energy of the corresponding 
integral circulant graphs would be most desirable. Of course, this should then be the basis 
for analyzing these graphs for hyperenergeticity. 
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